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A  MULTIPLE-ATTACK  EFFECTIVENESS  PREDICTION  ALGORITHM 


1.  INTRODUCTION 


When  launching  a  multiple  attack  against  an  evader,  there  are  several 
questions  that  must  be  considered.  Among  these  are:  When  is  the  best  time  to 
launch?  How  many  attackers  should  be  deployed?  What  are  the  optimal  attack 
trajectories?  To  begin  to  answer  these  questions,  a  real-time  computational 
procedure  is  needed  that  predicts  the  effectiveness  of  multiple  attackers 
against  a  prescribed  invader. 

One  measure  of  attack  effectiveness  is  given  by  the  radial  distance 
between  evader  and  attacker  at  a  specified  termination  time.  Varying  degrees 
of  damage  can  be  inflicted  by  simply  placing  the  attacker  in  proximity  to  the 
evader  at  that  time.  Accordingly,  range  from  attacker  to  evader  is  perhaps 
the  single  most  important  parameter  affecting  performance  in  any  particular 
situation.  Under  realistic  operating  conditions,  discrepancies  between  the 
actual  and  estimated  values  of  range  will  invariably  occur  because  of  attacker 
delivery  errors,  as  well  as  errors  inherent  in  localizing  the  evader.  To 
reliably  predict  attack  performance,  these  errors  must  be  accurately  modeled. 

Such  a  model  is  developed  analytically  in  this  report.  The  model  allows 
the  probability  of  evader  disablement  to  be  computed  as  a  function  of  attacker 
characteristics,  evader  characteristics,  and  associated  error  statistics.  The 
factors  that  significantly  contribute  to  overall  attack  effectiveness  are 
quantified.  Typical  inputs  include  (1)  evader  localization  error  statistics 
and  the  associated  position  and  velocity  estimates,  (2)  evader  evasion 
characteristics,  (3)  attacker  trajectory  error  statistics,  (4)  attacker 
destruction  characteristics,  (5)  mean  evader  position  after  alertment,  and  (6) 
an  offset  between  the  mean  evader  position  and  attacker  termination  point. 

The  model  is  subsequently  transformed  into  a  compact  computational  algorithm 
suitable  for  making  online  attack  performance  predictions. 

In  section  2,  the  problem  is  formulated  and  equations  for  the  probability 
of  evader  disablement  are  developed.  Section  3  provides  simplifying 
approximations  to  the  probability  distributions,  and  section  U  presents  a 
closed-form  expression  for  the  probability  of  evader  disablement.  In  section 
5,  the  results  are  modified  to  include  multiple  attacks  launched  against  a 
single  evader. 
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2.  PROBLEM  FORMULATION 


Consider  the  geometric  description  of  an  attacker  and  evader  encounter 
depicted  in  figure  1.  Initially,  the  attacker  is  launched  on  a  trajectory 
specified  by  the  vector  X^  *  (X*  Yw),  where  the  prime  symbol  denotes 
matrix  transposition.  However,  because  of  inherent  errors  in  the  delivery 
system,  the  attacker  will  not  necessarily  impact  at  its  intended  termination 
point.  To  accommodate  possible  discrepancies  between  these  two  points,  the 
actual  attacker  position  at  time  of  termination  is  denoted  by  X^  *  (Xw,  Yw). 

For  computational  convenience,  all  errors  associated  with  attacker  positioning 
are  tacitly  assumed  to  be  Gaussian-distributed  with  the  following  statistics: 

EtX*]  -  X*  ,  (la) 

E[(Xw  -  X^QC*  -  X*)*]  .  Pw,  (lb) 

where  E[»]  is  the  statistical  expectation  operator. 


•cow  ms  cvAotn  localization  errors 

AMO  EVASION  CHARACTERISTICS 


Figure  1.  Geometric  Description  of  Attacker/Evader  Encounter 


In  a  similar  manner,  the  actual  and  estimated  values  of  evader  position 
at  time  of  termination  are  defined  by  X^  *  (Xt,  Yt)  and  X^  =  (Xt,  Yt), 
respectively.  To  determine  the  statistical  characteristics  of  Xt,  observe  that 
(see  figure  2) 

Xt  -  Xto  +  Xta  ,  (2) 

where  Xto  represents  evader  position  at  alertment  of  an  attack.  Estimates 
of  the  vector  Xt0  are  typically  obtained  from  localization  and  motion 
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BEFORE  ALERTMEMT 


Figure  2.  Vector  Diagram  of  Relative  Evader  Position  at  Termination 


analysis  algorithms,  which  may  also  provide  statistical  descriptions  of  the 
associated  estimation  errors.  Since  these  errors  are  modeled  as  Gaussian 
random  variables  (reference  1),  Xt0  will  also  be  Gaussian-distributed,  with 
mean  and  covariance  given  by 


E[Xto] 

*  -to  » 

(3a) 

E[(Xto 

■  ?to^^to  ■  ?to^'^  a  ^to  * 

(3b) 

The  vector  X^a  defines  the  evader  position  at  time  of  termination 
relative  to  its  position  at  time  of  alertment.  Unfortunately,  errors 
associated  with  estimating  Xta  have  a  complicated  non-Gaussian  probability 
distribution  that  functionally  depends  on  the  evasion  characteristics.  To 
devise  tractable  approximations  for  this  distribution  requires  a  suitable 
geometric  model  of  the  evader  escape  region.  The  simplest  model  is  a  circle 
centered  at  Xt0  whose  radius  increases  with  time  in  proportion  to  the 
maximum  evader  speed.  Although  such  models  are  frequently  used  in  practical 
applications,  they  are  somewhat  unrealistic  because  they  implicitly  assume 
that  the  evader  has  a  zero  turning  radius.  Recently,  more  realistic 
techniques  of  modeling  evader  evasion  characteristics  have  been  developed 
(reference  2).  The  resulting  escape  regions,  shown  in  figure  3,  accurately 
reflect  constraints  on  the  evader  dynamics  by  utilizing  a  non-zero  turning 
radius  in  conjunction  with  a  finite  maximum  speed.  At  any  prescribed  time 
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Figure  3.  Growth  of  Evader  Escape  Region  with  Time 


after  alertment,  the  evader  must  reside  within  a  region  whose  boundary 
consists  of  an  involute  (dotted  curve  in  figure  3)  and  arcs  of  the  two 
minimum-turning-radius  circles. 

As  an  example,  the  darkest  shaded  area  in  figure  3  depicts  the  evader 
escape  region  at  1  minute  after  alertment.  In  this  illustration,  the  origin 
of  coordinates  is  defined  by  the  vector  Xt0  (the  assumed  evader  position  at 
alertment). 

Because  the  actual  escape  region  is  difficult  to  describe  analytically, 
it  will  be  approximated  by  a  circumscribing  circle  as  shown  in  figure  4. 
Although  this  approximation  introduces  errors  into  subsequent  probability 
cosiputations,  these  errors  decrease  with  time  because  the  actual  and 
approximate  escape  regions  eventually  coincide  (reference  2).  For  convenience 
of  notation,  R4  is  used  to  denote  the  radius  of  the  circumscribing  circle, 
and  Xta  defines  the  location  of  its  center  relative  to  Xt0. 

While  it  is  known  that  an  alerted  evader  must  necessarily  reside  within  a 
prescribed  escape  region,  more  specific  data  to  further  refine  the 
localisation  process  are  usually  unavailable.  Evasive  maneuvers  are  almost 
certain  to  occur  because  of  the  crisis  nature  of  the  situation  (i.e.,  the 
evader  knows  it  is  under  attack).  However,  the  absence  of  additional 
information  precludes  classifying,  in  general,  one  type  of  maneuver  as  more 
likely  than  another.  In  light  of  this  fact,  it  is  not  unreasonable  to  assume 
that  the  errors  associated  with  estimating  Xta  are  uniformly  distributed 
over  the  prescribed  escape  region.  Under  such  an  assumption,  it  can  be  shown 
that 
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Figure  4.  Circular  Approximation  of  Evader  Escape  Region 


E(Xta] 


5ta  • 


E[(Xta  "  XtaH*ta  -  *ta>']  ■  (V*)2  I  *  <*b) 

where  I  ia  the  two-dimensional  identity  matrix,  and  the  pertinent  statistical 
computations  are  performed  using  a  circular  approximation  as  depicted  in 
figure  4. 

To  determine  the  probability  of  evader  disablement  resulting  from  an 
attack,  it  is  necessary  to  mathematically  combine  the  statistical  information 
embodied  in  equations  (1)  through  (4)  with  known  attacker  damage  character¬ 
istics.  This  task  may  be  accomplished  via  application  of  the  law  of  total 
probability  (reference  3),  which  leads  to  the  expression 

ProblEI]  -  J  Prob[El|R]p(R)dR,  (5) 

Q 

where  R  *  |x|  is  the  range  from  evader  to  attacker  at  time  of  termination, 
and  Q  is  the  set  of  all  possible  values  of  R.  Here,  the  conditional 
probability  function,  Prob[El|R],  defines  the  probability  of  evader 
disablement  given  that  the  relative  range  from  attacker  to  evader  is  known  at 
time  of  termination.  It  is  this  term  that  incorporates  attacker  destruction 
characteristics  into  the  effectiveness  prediction  model.  The  remaining  term 
p(R)  represents  the  probability  density  of  R  and  is  used  to  account  for  evader 
localization  errors  and  evasion  characteristics,  as  well  as  attacker 
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trajectory  srrors.  To  see  this,  observe  from  figure  2  that  the  vector  X  can 
be  written  as 

X  =  Xto  +  Xta  -  X*  •  (6) 

Thus,  both  the  magnitude  and  direction  of  X  is  affected  by  evader  localization 
errors  via  Xto,  evasion  characteristics  via  Xta,  and  trajectory  errors  via  X*. 

To  compute  the  probability  of  disablement  online  and  in  real  time,  it  is 
first  necessary  to  explicitly  represent  p(R)  and  Prob[El|R]  as  integrable 
functions  of  R.  Unfortunately,  it  is  virtually  impossible  to  develop  exact 
mathematical  representations  of  these  two  functions  that  are  amenable  to 
closed-form  integration.  Nevertheless,  suitable  approximations  can  be 
employed,  allowing  equation  (5)  to  be  evaluated  without  recourse  to  numerical 
integration.  The  desired  approximations  for  p(R)  and  Prob[El|R]  are 
computed  in  the  next  section. 
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3.  APPROXIMATION  OF  THE  PROBABILITY  DISTRIBUTION 


To  mathematically  describe  the  probability  density,  p(R),  it  is 
convenient  to  utilize  the  concept  of  characteristic  functions  (reference  3). 

For  the  problem  under  consideration,  observe  from  equation  (6)  that  the  vector 
X  is  a  linear  function  Xt0,  Xta,  and  X*.  Since  the  errors  associated 
with  estimating  these  latter  three  vectors  are  uncorrelated,  it  follows 
(reference  3)  that  the  characteristic  function  of  p(X)  satifies  the  relation 

4>X(«)  -  *to<»><*ta<«>,M-«>  •  (7) 

where  ^ta^St^*  an<*  are  characteristic  functions  of  the  probability 

densities  p(Xt0),  p(Xta),  and  p(Xw),  respectively.  In  view  of  the  assumption 
that  Xt0  and  X*  are  Gaussian-distributed,  the  functions  <&to(«)  and  4>w(w) 
take  the  form  (reference  3) 

*to<«)  ■  exPtj“'Hto  ~  (l/2)«,Pt0«]  *  (8) 

$w(m)  *  exp(joi'Xw  -  (l/Zju'Pyu]  .  (9) 

As  discussed  previously,  the  vector  Xta  is  presumed  to  be  uniformly 
distributed  over  the  evader  region;  thus 
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and  L^P)  [•]  are  Laguerre  polynomials  (reference  4)  defined  by 


LiP)m 

m 


m,p 


0,1,2 


( 


Here,  n  appears  as  an  argument  in  the  characteristic  function  of  pn(^ta) 
to  denote  the  order  of  the  approximating  series. 

Since  equations  (7),  (8),  (9),  and  (13)  completely  specify  the 
characteristic  function  of  p(X),  this  probability  density  may  be  determined 
analytically  as  follows: 


p(X)  =  — 4>  (<o)  exp[-ju'X]dw  .  ( 

{It)1  X  "  ~ 

To  obtain  the  distribution  of  range  errors,  a  transformation  to  polar 
coordinates  is  made,  such  that 


p(R,9)  *  p(X) | | J | | 


where 


— i  *v(tt>  exp[-juR  cos  (0  -  <J>)]du  ,  ( 

(2ir)2  X 


*  tan~l  (ojy/ujj)  .  ( 

The  desired  probability  density  p(R)  follows  by  integrating  p(R,9)  over 
the  interval  0  <  0  <  2ir;  that  is, 

p(R)  -  oJ2irp(R,0)d0 

■  — ^ /°_jT  <>x(«)  [o^2^  exp[-jwR  cos  (0  -  <t>)]Jdu 

-  h  -O’  V^)R  V“R><*“  ’  ( 

where  Jg(o)R)  is  a  zeroth  order  Bessel  function  (reference  5). 

To  mathematically  describe  the  function  Prob[El|R]  requires  explicit 
knowledge  of  the  attacker  destruction  characteristics.  For  this  study, 
attackers  with  circular  interdiction  zones  are  used. 

The  center  of  the  zone  coincides  with  attacker  position  at  time  of 
termination.  For  computational  convenience,  R4  will  be  used  to  define  the 
radius  of  the  destruction  envelope. 
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Although  specific  numerical  data  are  available,  an  exact  mathematical 
description  of  P [ El | R]  does  not  exist.  However,  from  practical  considerations, 
it  is  apparent  that  this  function  will  decrease  monotonically  as  R  increases 
and  also  must  satisfy  the  inequality 


0  <  P[EI|R]  <  1  , 

for  0  £  R  ®.  These  characteristics  suggest  that  finite-sum  Gaussian 
distributions  again  be  used  as  a  basis  for  approximating  P[El|R],  i.e., 

P(EI|R]  «  P[EI | R,m] 

where 

.  m-1  k  „ 

P[El|R,m]  *  exp[-m(;  ]  l  rr  5 

k*0 

As  before,  m  determines  the  order  of  the  approximating  series,  and 

<;  -  R/Rd  • 


(21) 


(22) 


(23) 


(24) 


Graphs  of  equation  (23)  for  m  *  1,2,3,...,  10  and  m  =  ®  are  plotted  in 
figure  5.  As  expected,  these  functions  all  decrease  monotonically  as  R 
increases.  Note  also  that  the  limiting  conditional  probability  may  be 
expressed  in  the  form 


PIE1  |R,“]  *  bjS/*'1-  <25> 

Because  of  its  simplicity,  equation  (25)  is  frequently  used  in  practical 
applications  even  though  it  is  somewhat  artificial.  Indeed,  such  an 
approximation  implicitly  assumed  that  attacks  will  have  uniformly  devastating 
effects  over  the  entire  attacker  destruction  envelope.  In  light  of  this,  it 
is  perhaps  more  realistic,  and  certainly  more  economical,  to  employ  lower- 
order  approximations  when  possible.  Accordingly,  the  parameter  m  will  not  be 
specified  here  so  as  to  enhance  modeling  flexibility  and  realism. 


Figure  5.  Plots  of  the  Function  P[El|R,m] 
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4.  COMPUTATION  OF  THE  EVADER  DISABLEMENT  PROBABILITY 


The  probability  of  evader  disablement  resulting  from  an  attack  may  now  be 
determined  analytically  by  utilizing  approximations  developed  in  section  3  in 
conjunction  with  the  law  of  total  probability.  More  precisely,  substituting 
equations  (20),  (22),  (23),  and  (24)  into  equation  (5),  and  evaluating  the 
resultant  integral  over  Q,  leads  to  the  expression 


P[EI]  ■  $x(«)f(w)dM  , 


where 


«->  ■  r,  Ci 


[k“2]  '»p[|pd“2]  . 


Pd  35  « 


(26) 

(27) 

(28) 


(1) 


and  L'*'[»]  is  the  Laguerre  polynomial  defined  previously.  In  deriving 
equation  (26),  it  might  be  noted  that  integration  of  the  exact  probability 
distributions  over  Q  is  mathematically  equivalent  to  integrating  the 
approximate  distributions  over  a  semi-infinite  interval  0  <  R  <  ®. 


If  equations  (7),  (8),  (9),  (13),  and  (27)  are  combined,  the  integrand  in 
equation  (26)  may  be  rewritten  as 


-  si  Ci  [k*1]  [k*2]  -  kp“'  • 


(29) 


where 


■  X„  +  X„  -  X  , 
-to  -ta  -w  ’ 


>•[:] 

r?  (i,i)  p  d,2)i 
Lp  (1,2)  P  (2,2)J 


(30) 


Pto  +  Pw  +  (pe  +  pd)X  • 


(31) 


and  I  is  a  (2  x  2)  identity  matrix.  Utilizing  equation  (16),  in  conjunction 
with  equations  (26)  and  (29),  then  yields 


PlEI‘  '  *  Jo  llo  P«  • 


where 


k+1 


s+i 


L_  i=ll 

2w  k!  i! 


0 


2\k+i 


explju'X  -  Pu]dw 


(32) 


(33) 
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Th«  reaeining  tack  ic  to  determine  in  terns  of  the  pertinent 

Modeling  parameters.  This  nay  be  accomplished  by  employing  equation  (14)  and 
the  binoaial  expansion  theorua  to  explicitly  evaluate  the  integral  appearing 
in  equation  (33).  After  integrating  and  simplifying,  equation  (33)  takes  the 
fora 


♦  OJ  /k*i\  k+l 

^  U  )  Jo  9k*l-ll  VV  ’ 


9,u.x]  -  (?)  [-*xr*  \jfi  . 

where  X*  -  are  the  eigenvalues  of  P,  given  by 


I —  I  «  II  x ~ i  f - •  \-i 

2  J  VL  2 


H  .  P  (1,2) 


♦  ■  tan 


-if  -  10x11.  1 

Lxx  -  p(l*l)J 


Substituting  equation  (34)  into  equation  (32),  the  probability  of  evader 
disableaent  for  one  attacker  is  deterained  as  a  function  of  attacker 
destruction  characteristics,  attacker  trajectory  coordinates  and  associated 
error  statistics,  evader  localisation  solution  and  error  statistics,  or  evader 
evasion  characteristics. 


5.  MULTIPLE  ATTACK  PREDICTION 


The  previous  analysis  determined  the  probability  of  evader  disablement 
for  a  single  attack  (equation  (32)).  In  this  section,  the  combined 
probability  of  evader  disablement  is  developed  when  more  than  one  attacker  is 
launched  at  a  single  evader.  To  reduce  the  complexity  of  the  algorithm,  only 
exponential  destruction  characteristics  (m  ■  1)  are  employed,  and  the 
termination  times  for  all  attackers  are  assumed  to  be  identical  as  well  as  the 
destruction  radii.  In  addition,  trajectory  errors  are  assumed  to  be  zero. 

From  equation  (23),  the  conditional  disablement  probability  distribution  for 
exponential  destruct  characteristics  is  given  by 


P (BI | R , 1 ]  -  exp[-$2]  , 

(39a) 

i  -  R|Rd  • 

(39b) 

R  *  l?to  +  Hta  “  5wl  * 

(39c) 

where  R  is  the  range  from  attacker  to  evader  at  termination. 


To  determine  the  probability  of  evader  disablement  for  multiple  attacks, 
modifications  to  the  problem  formulation  in  the  previous  sections  are 
necessary.  Originally,  evader  disablement  probability  was  found  by 
integrating  equation  (5)  with  respect  to  R  (range  from  evader  to  attacker  at 
termination).  However,  for  multiple  attackers,  there  is  a  termination  range 
for  each  attacker.  To  alleviate  this  problem,  trajectory  errors  are  assumed 
to  be  zero.  Then,  equations  (6),  (7),  and  (17)  reduce  to 


a  « 


*to  -ta 


(AO) 


P(X)  -  — r  V«)  «xPU»'XJd“  » 

~  (2w)2  X  * 

♦x(m)  ■  ♦tote^tate)  • 

For  multiple  attacks,  there  is  a  conditional  probability  distribution  for 
each  attacker  that  takes  the  form 

PtBllRj  -  exp[-cj]  ,  (A3a) 

P[EI|R2]  ■  exp[-C2]  ,  (A3b) 

PlEllRj,]  -  exp[-*2]  ,  (A3c) 


where  | and  R^,  R£,  . ...  Rjj  are  the  ranges  from  the  evader 

to  each  attacker.  If  there  is  no  overlap  of  the  destruct  regions,  the  total 

probability  for  evader  disablement  is  the  sum  of  the  disablement  probabilities 


from  the  individual  attackers  given  in  equation  (32)  and  can  be  written  as 


Pt(EI]  =  Pi(EI]  +  P2[EI]  +  ...  +  Pn[EI]  .  (4 

When  there  ia  overlap  of  the  regions,  which  realistically  occurs  due  to 
the  exponential  characteristics  of  the  attackers,  the  individual  conditional 
probabilities  must  be  combined.  A  formulation  well  suited  to  this  situation 
is  given  by 

N  9 

PtEI |R1,  R2 . Rfl]  -  1  -  n  [1  -  exp(-q)]  .  (4 

'  '  i*l 

Substituting  equation  (43)  into  (45)  and  expanding  yields 

N  ,  N=1  N-i  ,  , 

P(El|Rlt  R2,  ....  Rn)  »  l  exp[-qi  -  l  l  exp[-(q  +  .]] 

c  i»l  i»l  jal  ^ 

N-2  N-i  N-i-j  ,  -  , 

i  - 

+  (-UN_l  expMsJ  +  Cj  ♦  ...  t  5^)]  .  (4 

Recalling  that  R  ■  |x|  and  employing  equations  (6)  and  (24),  each  exponential 
for  the  attack  destruction  distribution  can  be  written  in  the  form 


...2  .2  ..  /pst  -m1*!2 

expt-Cfj  *  |  +  ...  )]  =  «xp(-  - -j -  ♦  - -J - 

1  Rj  R . 


where 


?t  *  Xto  +  ?ta» 


and  Xy(*)  is  the  position  of  attacker  (*)  at  termination.  Now  equation 
(46)  can  be  written  as 


,  »  r  ist-x^uji2 

PtEl|Xt,  X/l),  ^(2),  ....  X/N)]  *  l  exp  -  — - == - 

i=l  Rj 

a 

M;lN;1  /P*t is.  - 

-  I  l  exp/-  - = -  +  - 5 -  I... 


and  the  probability  is  given  by 


plEIl  -  _a,;a0_aB/a0  <J>x(«)f(«)do»  ,  (49) 


?(«)  -  — i-r  J* exp[-j «’XJ  PlEljX.,  X  (1),  ....  X  (N)JdX.  .  (50) 

(2w)2  1  t  -w  ^  t 

As  can  be  seen  from  equation  (48) ,  the  first  summation  weights  the 
destruction  capabilities  for  each  attacker  individually  in  the  exponential; 
the  double  summation  averages  the  destruction  capabilities  in  combinations  of 
two  in  the  exponential;  and  the  triple  summation  averages  the  destruction 
characteristics  in  combinations  of  three.  This  process  continues  until  the 
last  term  averages  all  the  destruction  characteristics  inside  the  exponential. 
Equation  (48)  can  be  written  in  the  form 


p[El|Xt,  X^l),  ....  X^OO]  *  l  (-1)1=1S.  ,  (51) 

ial 


where  S.  represents  the  first  summation  which  contains  (.)  exponential 

‘  N1 

functions,  represents  the  double  summation  and  has  (^)  exponential 

functions,  etc.  In  general  for  N  individual  weapons,  S0  represents  the 

N  * 

2th  summation  which  has  (^)  exponential  functions  where  2  *  1,  2,  N. 
Substituting  equation  (51)  into  (50)  results  in  N  integrals  of  the  form 


t.(«)  -  -  l~z  //  exp[-ju’X  ]  SjjdX  ,  (52) 

1  (2ir)2  2  fc 

where 

N  N  1 

?(«)  -  i  (-ir  lv«)  .  (53) 

2«1  * 

The  solution  for  t^Cm)  through  ^(u)  is  the  same  except  for  the  number  of 
attackers  contained  in  the  exponential.  The  integral  is  solved  for  %(u) 
in  the  remainder  of  this  section  and  then  the  solutions  for  ^(u)  are 
given.  Substituting  for  Sjf  in  equation  (52),  ij^w)  is  written  as 


17 


Let 


Rj 

-  d 

-  ■  r 


«  • 


■•■0 


X  (i) 


0 


then  equation  (54)  becomes 
1 


y«) 


».,* -r-r “p(- Y't  I-* Xi ' ^]) 


•  "p(- Y'%  Uyi'%]) 
(-  i*  LL  <*‘  *  r‘J)- 


dxdy 


exp 


From  reference  6,  the  integral  becomes 
«2 


(55) 


V*>  *  ss  exp 


~1  “ 


L^d 


1  N  2  2 
-  h  x.  <*i +  ❖ 


R“  i-1 

u 


•  exp 


where 


^  N  "  2^  2N  -J  ’ 


(56) 


(57) 


Solutions  for  ?,  («)  through  .  (u)  can  be  found  in  a  similar 
l  N-i  ^ 

manner.  For  t^(u),  the  solution  can  be  reduced  to  (^)  exponential 
functions  of  the  fora 
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I 

« 

. 

« 

i 

I 

I 

t 

1 


(58) 


rj  n  r  i  Rd  1 

V»>  -  expj-js.*^  -  2«*f  «J  . 

|| 

Similarly,  ?2(«)  has  (’)  exponential  functions  and  is  given  by 


*2<«> 


2(4*) 


N-l  N-i  /[X  (i)  + 

l  l  exP  (  - - 

i»l  j-1  \ 

2  2  \ 

IZiHin) 

*  / 


X  (i+j)]'  [X  (i)  +  X  (i+j)] 


2R 


2  2 
Xi  *  *1+1 


where 


(59) 


(60) 


Following  this  procedure,  the  remaining  solutions  for  tj/w)  can  be 
found  and  substituted  into  equation  (53)  to  obtain  t(»).  Substituting  for 
t(m)  into  equation  (49),  the  probability  of  evader  disablement  is  given  by 


PfEIJ  -  I  (-1)1"1  J-"  /"  Vu)t,(«)d«  • 

i-l  *  1 

which  can  be  written  as 

N  i-l 

P[EI]  •  Z  (-1)  P.IEIJ  , 
i-l  1 


where 


PjlEI]  - .jJ* _J>m  ♦J{(m)ti(«)d»  . 


(61) 


(62) 


(63) 


Once  again,  the  solution  [El ]  through  Pj|[EI]  is  the  same  except  for 
the  number  of  weapons  contained  in  the  exponential.  A  solution  for  P^lEI] 
is  given,  and  then  the  remaining  solution  can  be  found.  Substituting  equations 
(8),  (13),  (15),  (16),  (42),  and  (56)  into  equation  (63),  PN[EI]  can  be 
written  in  the  form 
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K .  n“l  /  \ 

PN[EI1  *  N(2n)  jJq  (k+l)* 


where 


*  ■  i*  • 

°k*  ■  !*£■  -.O'  (r)  -  kpx^ld“  • 


"  -to  +  -ta  -  N  * 


PXN  '  Pto  *  V  *  25  1  ' 

The  solution  for  Qrn  ia  aimilar  to  the  solution  for 
(equation  (34))  and  takes  the  form 

exp[-  |(IL  ♦  £L)]  k 

^  UO  0k-!l 

where  ©£•*•]  is  given  by  equation  (35),  Xyy  and  Xyji  are  the  eigenvalues 
of  ?xn  as  in  equation  (36),  and  the  ■  (lxN>  '  ia  given  by 

f  “*  *m  ,ln  M  r;  i 


~7TmT  yXYN  ^ 

sin  ^  cos  ♦jj 

Ly» 


♦„  .  Un' 


■i  r  v1'21  i 
l>  -  p*(i-i,J 


Similarly,  P^lEI]  can  be  written  in  the  form 

'i'11'  ■  rfo  A  (k"i)  "I  “«  • 


l 


1 


mmw*  mwn  v.w*  n."s» 


°kl 


hi 


N 

l 

i-1 


expt-  |(£xi+ 

VXXi\iT 


*pL  9a^Xi* 

P-0 


\(i'  9k-SI5Yi’ 


\i] 


io  ♦  *«  -  V0*  1 " l-  *• 


P,  ,H.? 
to  e  2 


'•♦I 


N 


(72) 


(73) 


(74) 


[cos  sin 


lli 


~7Tn  ~^r 

sin  cos 

vxyi 


-li  * 


(75) 


tan 


-if  V1-21  1 

|>  -  fi(lil)J  ’ 


(76) 


and  ©[•*•!  is  given  by  equation  (35)  and  the  eigenvalues,  XXi  and  Xyi* 
are  given  by  equation  (36).  The  equations  for  P2(EI)  •••  ?N— 1 CEI )  can  be 
fomulated  in  a  similar  manner  and  then  substituted  into  equation  (62)  to  obtain 
the  total  probability  of  evader  disablement  for  N  attackers. 

The  previous  section  has  shown  that  to  determine  the  probability  of 
evader  disablement  for  multiple  attackers,  the  conditional  evader  disablement 
probability  breaks  down  to  sums  of  exponentials  as  is  shown  in  equation  (48). 

For  N  attackers,  each  summation  averages  the  destruction  characteristics  for 
each  attacker  individually  and  in  combinations  of  2,  3,  4,  ...,  N.  The  end 
result  is  the  total  probability  of  evader  disablement  is  determined  by 
integrating  each  exponential  in  all  the  sums  over  the  entire  uncertainty/ 
evasion  region  (Xto  +  Xta).  The  attacker  trajectory  errors  can  now  be 
included  by  determining  an  average  trajectory  error  in  each  summation.  The 
necessary  modifications  to  the  algorithm  are  presented  for  the  Nth  summation 
in  equation  (51).  Modification  to  the  other  summations  follow  in  a  similar 
fashion.  For  the  Nth  summation,  the  average  trajectory  error  can  be  written 
as 


Ks  “  N(^wl  +  ^w2  + 


+  X  M) 
-wN 


(77) 


where  X**  are  now  Gaussian-distributed  random  variables.  Since  the 
trajectory  errors  are  assumed  to  be  independent  Gaussian  random  variables, 
Xy||  is  also  Gaussian  with  the  following  statistics: 


ElLJ 


*WN 


1  N  - 

N  1  *wi 
"  i-1 


(78a) 


21 


and  the  characteristic  function  for  becomes 


♦wN(»>  -  - 1  “,pwN“i 


Now  equations  (40),  (41),  and  (42)  can  be  written  in  the  form 


s»  -  s,  -  s* 


P<V  -  ^2  -O*  ^(a)  «pHa’<xt  -  !„,)]*. , 


#^(«)  ■  ♦to(«)*ta(«)^wN(-<!i.)  * 


and  equation  (54)  is  now  given  by 


’  ,K<">  "  -  S*>1 


/PSt 

•  expi-  - 


-5wil 


ISt  -  *w2l 


♦  .  •  •  ♦ 


l?t  -  ^wNl 


dXfc  . 


With  the  addition  of  the  X„n  term  in  equation  (83),  and  recognizing 
that  3n/N  *  XvN>  the  solution  for  fjj(«)  reduces  to 


R?  L'  L  ,  N  ,  9 

M«)  *  rr™  exp  — j - «  l  <xi  + 

N  4irN  NR?  R?  i-1  1  1 

u  a 


?)  •  eXP[4  tt'ow  ttl 


The  next  equations  that  are  modified  due  to  the  addition  of  the 
trajectory  errors  are  equations  (66)  through  (70),  which  become 


^N  "  —to  +  -t.  -  *wN  * 


PXN  "  Pto  +  PeJ  +  2N  1  +  PwN  ’ 


°kN  " 


exp[-  |(fiL  ♦  k 

uo  01^xn’*xn1  ^-h^yn’^yn* 


wwv.v\n  v.wv VV-VLV.V 


-  •  »  ■  >  »>  T 


wfwwjuwuiwuwiwwiflumwuiwHHHwnmwnwwmuwwr 


wwwwmiwwwunmmiiinyiwwwiH 
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cos  $N  sin  $N  H 


i  in  $, 


N 


cos 


* 


N 


L  ^  J 


k 


-1 


p^a.2) 


*XN  "  ?XN(1,1) 


(88) 


(89) 


and  Xxn*  ^YN  are  the  eigenvalues  for  ?XN‘ 

Following  this  procedure,  the  other  equations  can  be  modified  and  the 
total  evader  disablement  probability  can  be  determined. 
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6.  CONCLUSIONS 


This  report  has  presented  a  real-time  computational  procedure  for 
predicting  the  effectiveness  of  an  attacker  against  a  single  evader.  A  j 

functional  relationship  was  developed  that  expresses  the  probability  of  evader 
disablement  in  terms  of  attacker  trajectory  and  associated  error  statistics, 
attacker  destruction  characteristics,  and  evader  localization  error  statistics 
and  evasion  characteristics.  Furthermore,  the  computational  algorithm  is 
sufficiently  compact  to  allow  for  continuous  updating  of  attack  performance 
predictions  as  tactical  conditions  change.  j 

Finite-sum  Gaussian  distributions  were  used  to  model  the  errors 
associated  with  estimating  evader  position  after  alertment  and  the  attacker 
destruction  characteristics.  By  selecting  the  appropriate  index  (n,m)  for 
each  of  these  finite-sum  series,  uniform  to  Gaussian  density  functions  can  be  j 

obtained  to  model  the  pertinent  distributions.  These  two  indices  are  left  as  | 

inputs  to  the  effectiveness  algorithm  so  as  to  enhance  modeling  flexibility 
and  realism. 

Finally,  the  probability  of  evader  disablement  was  extended  to  include  ] 

multiple  attacks.  To  reduce  the  complexity  of  the  algorithm,  only  exponential  i 

destruction  characteristics  (m  *  1)  were  modeled.  In  addition,  it  was  assumed  | 

that  the  termination  time  and  destruction  radii  for  all  the  attackers  were  the 
same.  The  trajectory  errors  associated  with  each  attacker  was  also  averaged 
into  a  single  error.  Total  mission  effectiveness  was  determined  by  employing  ! 

the  "additive"  probability  law.  An  algorithm  that  predicts  total  disablement 
probability  for  N  attackers  was  subsequently  presented. 

i 
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